Temperature-induced gap formation in dynamic correlation functions of the 

one-dimensional Kondo insulator 
— Finite-temperature density-matrix renormalizat ion- group study — 
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Combination of the finite-temperature density-matrix renormalization-group and the maximum 
entropy presents a new method to calculate dynamic quantities of one-dimensional many-body sys- 
tems. In the present paper, density of states, local dynamic spin and charge correlation functions 
of the one-dimensional Kondo insulator are discussed. Excitation gaps open with decreasing tem- 
perature and the gaps take different values depending on channels. The excitation spectra change 
qualitatively around the characteristic temperature corresponding to the spin gap which is the lowest 
energy scale of the Kondo insulator. 
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It is not rare that experiments on dynamic 
quantities such as photoemission spectrum, neu- 
tron inelastic scattering and optical conductiv- 
ity play a key role to understand the physics of 
strongly correlated electron systems. From theo- 
retical point of view, simple mean-field-type the- 
ories often fail to capture the essential physics of 
the strongly correlated systems and some unbi- 
ased method is required. This is the reason why 
many interacting systems have been studied in- 
tensively by various numerical methods includ- 
ing the quantum Monte Carlo simulations and 
the numerical exact diagonalization. 

The density-matrix renormalization-group 
(DMRG) method is one of the most powerful nu- 
merical algorithms recently developed fl]]. One 
of the authors and Wang and Xiang indepen- 
dently showed that the thermodynamic proper- 
ties of the ID XXZ-Heisenberg model are ob- 
tained by applying the DMRG method to the 
quantum transfer matrix, which may be called 
finite-T DMRG §1. 

Numerical study of dynamic quantities is more 
elaborated compared with static ones. The ap- 
plication of the maximum entropy (ME) method 
opened a way to obtain dynamic quantities reli- 
ably based on quantum Monte Carlo simulations 
1^, after various attempts had been made [§,0. 
However the problem to treat big enough sys- 
tems at low enough temperatures still remains. 
It is needless to say that the quantum Monte 



Carlo simulations suffer from the serious prob- 
lem of negative sign for most of fermionic sys- 
tems and quantum spin systems with frustra- 
tion. Calculation of dynamic quantities by the 
Lanczos algorithm of the exact diagonalization 
is based on the continued fraction [0. For this 
method the limitation of the system size is so 
severe that it is almost impossible to discuss ex- 
citation spectra in the thermodynamic limit. 

In this Letter, we show that the combination of 
the DMRG and the ME presents a new powerful 
method to calculate dynamic quantities of one- 
dimensional (ID) many-body systems at finite 
temperatures. Correlation functions along the 
imaginary time direction can be calculated by 
the finite-T DMRG method. After Fourier trans- 
formation, analytic continuation from the imag- 
inary frequencies to the real frequencies may be 
done through the Fade approximation p| . How- 
ever it has turned out that this procedure is of- 
ten unstable at low temperatures and we have 
found that the most efficient and reliable way to 
obtain the spectra is the ME method. The com- 
bination of the finite-T DMRG and the ME is 
free from statistical errors and the negative-sign 
problem. Since we deal with the transfer matrix 
we can discuss the thermodynamic limit directly. 
However, additional calculations are necessary to 
obtain nonlocal correlation functions and in the 
present paper we restrict ourselves to the local 
quantities. 
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As the first apphcation of the method we dis- 
cuss in this paper dynamic quantities of the ID 
Kondo insulators based on the Kondo lattice 
(KL) model. The KL is the simplest theoret- 
ical model for heavy fermion systems. By in- 
tensive studies in the last ten years, the ground 
state phase diagram of the ID KL model has 
been completed 0. At half- filling the ground 
state is always paramagnetic and insulating and 
therefore this phase is considered theoretical 
prototype of Kondo insulators in ID. 

The effect of the strong correlation in the 
ground state of the Kondo insulator appears in 
the difference between the spin gap Ag and the 
charge gap Ac p|,|lO|- It was shown by quantum 



Monte Carlo simulations and recently by the 



finite-T DMRG 121 that the difference between 



As and Ac are reflected in the temperature de- 
pendence of spin and charge susceptibilities and 
speciflc heat. 

In this Letter we discuss temperature depen- 
dence of the density of states, the dynamic spin 
and charge correlation functions. It is shown 
that the gap in the density of states develops be- 
low the smallest characteristic temperature cor- 
responding to the spin gap rather than the quasi- 
particle gap. Of course, the edge of the spectrum 
itself is given by the quasiparticle gap |TB|,|T^. It 
is a general feature of strongly correlated sys- 
tems that a change in the lowest energy scale 
can trigger a strong modiflcation of an excita- 
tion spectrum. Similarly, the charge excitation 
spectrum changes qualitatively at the character- 
istic temperature. Another consequence of the 
strong correlation is the two-body spectra of the 
dynamic spin and charge correlations are quite 
different from the convolution of the density of 
states. 

The Hamiltonian we discuss is the usual ID 
KL model deflned as 

H = -t E(4c.+i. + H.c.) + JY. S,i ■ Su, (1) 

is i 

where the operator Cis (cj^) annihilates (cre- 
ates) a conduction electron at site i with spin 

S ( = T,i)- Sci = J2a,l3cla{(^)apCii3/2 is the Spiu 

operator for conduction electrons and Su is the 
/-spin operator. The hopping matrix element 



is given by t, and J is the antiferromagnetic ex- 
change coupling between the conduction electron 
spins and the localized spins. The density of con- 
duction electrons is unity at half-flUing. In the 
present study, we treat the following dynamic 
quantities: the density of states p{uj), the local 
dynamic charge correlation function A^(ti;), the 
local dynamic spin correlation function of con- 
duction electrons Sduj) and that of the /-spin 

Here we outline the algorithm of the flnite- 
T DMRG method. In this method we use the 
quantum transfer matrix deflned by 



Tn{M) 



-Ph2n-l,2n/M -l3h2„,2n + l/M]M 



(2) 



where M is the Trotter number. In this for- 
malism the Hamiltonian H is decomposed into 

the two parts iJodd = En=l /i2n-l,2n and ifeven = 
L /2 

Yl,n=\ h2n,2n+i SO that wc cau evaluate the ma- 
trix elements of the exponential function. Since 
the partition function is given by the trace of 
the product of 7^(M), the partition function in 
the thermodynamic limit is determined from the 
maximum eigenvalue of 7^(M). In order to ob- 
tain the maximum eigenvalue and its eigenvec- 
tors for a large M, we iteratively increase M 
as follows. Note here that the right eigenvec- 
tor and the left one are different since 
the transfer matrix is not symmetric. We flrst 
diagonalize 7^(M) with a small M and obtain 
the eigenvectors corresponding to the maximum 
eigenvalue. These eigenvectors are used to con- 
struct the generalized density matrix. We di- 
agonalize the density matrix and use important 
eigenvectors which have large eigenvalues. Us- 
ing these eigenvectors as new bases we increase 
M within a flxed number of bases. Repeating 
the above procedures we obtain the maximum 
eigenvalue and its eigenvectors of %i{M) for a 
sufficiently large M at a given temperature T. 

The dynamic quantities are obtained from the 
correlation functions in the (3 direction. To get 
sufficient accuracy, it is necessary to use the 
finite-system algorithm of the DMRG By us- 
ing and obtained after the convergence 
of the finite-system algorithm, the spin correla- 
tion function, for example, is evaluated as 
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(3) 



-tH 



where ^^^(r) = e^^ Sffi 
tions will be published elsewhere 



Details of calcula- 
131 



The fol- 



lowing calculations are performed by the finite- 
system algorithm keeping 40 states per block. 
Truncation errors in the finite-T DMRG calcu- 
lation are typically 10~^ and 10~^ at the low- 
est temperature. The Trotter number used is 
M = 60, except for M = 80 at the lowest tem- 
perature. In this paper, we take t as units of 
energy and J is set to be 1.6. 

To obtain spectra of the dynamic quantities we 
solve the following type integral equations; 



dujSr(u!)e 



(4) 



We use the usual ME method with the clas- 
sic ME criterion and the flat default model 
PJT6[|. However comments are in order on how 
to rewrite the integral kernel and how to set sta- 
tistical errors for the ME method. 
Since there is a relation 



(5) 



with X{uj) being N{uj) or 5'c(f)(ci;), we use the 
kernel defined by 



cosh(| — t)uj 
coshf^ ' 



and rewrite the integral equation (^) as 

roo _ 

xK^j)= dujSf{uj)K{Tj,uj), 



(6) 



(7) 



1 _L Q-Pul 

Sfiuj) = —^S,{u). (8) 

As a consequence of using the kernel of @ the re- 
lation @ is satisfied with high accuracy; relative 
errors are much less than 1.0 x 10"^. It is possible 
to use the same kernel also for p{u) by using the 
relation p{uj) = p{—uj). Concerning the statis- 
tical errors aj for the imaginary-time data such 

as Xf(rj), we use aj = 1.0 x lO'^^-^Xf (T?) since 
truncation errors are the only source of errors in 



the finite-T DMRG method. We have confirmed 
that the spectra are not so sensitive for the de- 
fault models, the ME criteria, and the choice of 
a/s. 

Figure |1] shows the density of states p(u;) at 
various temperatures. At the lowest tempera- 
ture T = 0.1 shown in the inset, the gap com- 
pletely opens around the chemical potential u = 
0. The edge of the spectrum is consistent with 
the quasiparticle excitation gap Aqp = 0.7 ob- 
tained by the zero-temperature DMRG method 
|]T^,|T^ . The sharp peaks at the edges of the spec- 
trum indicate the formation of the heavy quasi- 
particle bands at low temperatures. The spec- 
trum has intensities also in the region of high fre- 
quencies, higher than the edge of the free conduc- 
tion electron band: u = 2. The high-frequency 
tail originates from the multiple spin excitations 
accompanied with the quasiparticle excitation. 
There is a dip between the sharp peak at the 
edge and the broad maximum at high frequen- 
cies. The dip structure indicates that spectral 
weight in the region is transfered, which may be 
similar in nature to the Fano antiresonance effect 
11. 
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FIG. 1. Spectra of single-particle excitations p{ui) at vari- 
ous temperatures. The inset shows the spectrum at the lowest 
temperature, T — 0.1. 



Let US discuss temperature dependence of 
p{uj). The most remarkable feature is that the 
quasiparticle gap is rapidly filled in with increas- 
ing temperature. At lower temperatures than 
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the quasiparticle gap, the gap structure aheady 
disappears and it becomes a pseudo-gap. At 
the same time the sharp quasiparticle peaks lose 
their intensities. The intensities are transfered to 
higher frequencies and the dip structure is also 
smeared out. At high temperatures the spec- 
trum shows only broad peaks which shift gradu- 
ally to ~ ±2. The broad peaks correspond to 
the band edges of the free conduction electrons. 
In the ID system, there are strong van Hove sin- 
gularities at the band edges uj = ±2. This be- 
havior is reasonable since the conduction elec- 
trons are gradually decoupled from the / spins 
as temperature is increased. 

Now we turn to the local dynamic charge cor- 
relation function N{(jj). A clear gap is developed 
in N{uj) at the lowest temperature. Fig. The 
spectrum of N{uj) at T = 0.1 has the edge at 
a; ~ Ac = 1.4, which is twice of Aqp = 0.7. 
There is a peak just above the gap. The peak 
structure corresponds to the charge excitations 
between the two sharp quasiparticle peaks of 
p{uj) shown in Fig. ^ The peak in N{uj) quickly 
changes into a shoulder as the temperature is 
increased. 
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FIG. 2. Spectra of charge excitations N{uj) at various 
temperatures. 



It is interesting to note that the intensity of 
N{u) is small above u; ^ 4, while the density of 
states p{uj) extends to much higher frequencies 
than uj = 2. This indicates that the structure of 
N{uj) cannot be understood by the simple pic- 



ture of renormalized band corresponding to p{uj). 
In contrast to the low-energy excitations, the 
high-frequency excitations around a; ~ 4 may 
be understood in terms of the excitations of the 
nearly free conduction electron band. The large 
intensity of the peak is due to the van Hove sin- 
gularities at the band edges of the ID system. 

At finite but low temperatures, there is a peak 
around ~ 0, which comes from the sharp 
peaks in p{uj). Quasiparticles and holes ther- 
mally populated at low temperatures contribute 
to the charge excitations within each peak. The 
integral intensity of the peak around ^ grows 
up to T ^ 0.2 as the temperature is increased. 
As the temperature increases further, the charge 
excitation gap in N{u}) is rapidly filled in. At 
T = 0.3 which is much smaller than the charge 
gap Ac = 1.4, one can no longer find any gap 
structure, similarly to the temperature depen- 
dence of p(uj). 

Next we show spin excitation spectra in Fig. 
Ij. The spectrum of Sc{uj) or Sf^uj) at the lowest 
temperature shows a clear gap and the gap edge 
is lower than that of N{uj) as expected. The gap 
edge is consistent with the spin gap Ag = 0.4 
obtained by the zero-temperature DMRG. 

At low temperatures T = 0.1,0.14, and 0.2, 
both Sc{uj) and Sfi^u) have two peaks. The 
peak at lower frequencies contain the /-spin ex- 
citations as the dominant component, while the 
main component of the peak at higher frequen- 
cies is the spin excitations of quasiparticles. In 
the high-frequency region of S'c(w), there is a 
shoulder structure at the lowest temperature, 
T = 0.1. This may be related with the fact that 
the lowest spin excitation has the momentum 
TT and therefore a significant part of the spec- 
tral weight of quasiparticle excitations around 
u; ~ 4 is transfered to the low-energy excitations 



T|,|T9| 



As the temperature increases, the higher- 
frequency peak of S{{u) and the lower- frequency 
peak of S'c(ci;) lose their intensities. Furthermore 
overall structures of N^u) and Sduj) are simi- 
lar at high temperatures and both have a 
peak structure around ~ 4. It shows that the 
main contribution to the peak is from the nearly 
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free conduction electrons. These results indicate 
that the conduction electrons and the /-spins are 
decoupled as the temperature is increased. 
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FIG. 3. Spectra of (a) Sc{i^) and (b) Sf{ui) at various 
temperatures; T = 0.10,0.20,0.30 and 0.40. The inset of 
(b) shows the comparison of two spectra, Sc{i-l>) and Si{uj) at 
T = 0.14. 



pie renormalized-band picture. One remarkable 
example is that the spectrum of A^(ci;) is differ- 
ent from the simple convolution of p{lj) partic- 
ularly at high frequencies. Another example in 
the ground state is the difference between the 
spin gap and the charge gap, which have been 
discussed previously 0. 

The temperature at which the quasiparticle 
gap and the charge gap disappear or the struc- 
ture of p{uj) and N{uj) vary drastically is T ~ 
0.3. This temperature corresponds to the spin 
gap Ag. In the present system, the singlet bind- 
ing between the conduction electrons and the /- 
spins produces the spin gap. As temperature 
is increased up to the order of Ag, the whole 
electronic states of the system are reconstructed. 
Therefore the smallest energy scale Ag influences 
also the temperature dependence of charge ex- 
citations. This is another new feature of the 
Kondo insulators observed for the first time by 
the present method. 
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In conclusion, we have calculated dynamic cor- 
relation functions of the ID KL model using the 
finite-T DMRG method and the ME method. 
At low temperatures the density of states p{uj) 
has a quasiparticle gap and the sharp quasipar- 
ticle peaks at the gap edges. The gap is rapidly 
filled in and the sharp quasiparticle peaks are 
broadened as the temperature is increased. The 
temperature dependence of p{uj) shows clearly 
the many-body origin of the gap formation. The 
spectrum of A^(ti;) has a charge gap at low tem- 
peratures, and the charge gap also disappears 
as the temperature is increased. It is not possi- 
ble to represent the many-body effects by a sim- 
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